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Abstract 

We give an algorithm to compute all the local peaks in the fc-level of an 
arrangement of n lines in 0(n log n) + 0{{kn)^^''') time. We can also find 
r largest peaks in 0(n log^ n) + 0((Tn)^^^) time. Moreover, we consider 
the longest edge in a parametric minimum spanning tree (in other words, 
a bottleneck edge for connectivity), and give an algorithm to compute 
the parameter value (within a given interval) maximizing/minimizing the 
length of the longest edge in MST. The time complexity is d{n^^''k^^'' + 



1 Introduction 

The A:-level of an arrangement of lines is one of popular geometric objects in 
computational geometry |12, 29|. The /c- level is the union of fc-th lowest (closed) 



line-segments of the arrangement, and it can be considered as the trajectory of 
the fc-th smallest element in a set of n data each of which depends on a parameter 
X linearly. Thus, the fc-level is a special case of the locus of the largest element 
of the minimum base of a parametric matroid ^ 19, ^ with one 



parameter; in precise, the fc-level is the locus of the maximum element in the 
minimum base of a parametric uniform matroid of rank fc. From a different 
aspect, the fc-level is a dual concept of the fc-set fl^ , and the complexity gk{n) 
of the fc-level of an arrangement of n lines is asymptotically same as the number 
of different fc-sets in a set of n points in a plane. 

Lovasz first gave a nontrivial 0{n^^'^) upper bound for 5[„/2j(f^), and 
also introduced Straus's fl(nlogn) lower bound. The current best upper and 
lower bounds for gk{n) are 0(fc^/^n) and 2^^'\/'°sfe)^ respectively. The 
upper bound holds for any parametric matroid (with a linear parameter) of rank 
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k in n elements 13 . Moreover, Cole et al. gave an 0{gk{n) log^ n) js] time al- 
gorithm to compute the fc-level; this time complexity has been improved due to 
improvement of algoirthms for the dynamic convex hull computation employed 
as a key subroutine ^. The current best time complexity (randomized ex- 
pected time) is 0{gk{n)a{gk{n))logn) time given by Har-Peled where a{) 
is the inverse Ackermann function. Thus, the time complexity is very close to 
0{gk{n) logn). 

However, we often need a compact "outline" of a trajectory of a parametric 
problem by using a small number of characteristic points on it. Such an out- 
line, generally speaking, will be useful as a compact data to control parametric 
problems, and possibly utilized in designing kinetic data structures jj, ^. Local 
peaks in the trajectory are considered to be natural characteristic points. A key 
observation to investigate the /c-level is that it has at most 2k — 1 local peaks (at 
most k maximal peaks and at most k—1 minimal peaks) ||, |l3|. One interesting 
problem is to compute all the local peaks efficiently. This enables us to give a 
decomposition of the fc-level into monotone chains, and hence create an outline 
with a size 0{k) of the fc-level. 

We give an algorithm to compute all the local peaks in the fc-level of an 
arrangement of n lines in 0{nlogn) + 0{{krif'^^) time, where O is the big-0 
notation ignoring polylogarithmic factors. The current estimate for the poly- 
logarithmic factor of the second term is less than log^ n; however we do not give 
it explicitly in this paper, since it is probably loose and will confuse readers. 
The time complexity is better than 0{gk(n) logn) for some restricted range of 
k even if the current lower bound of g kin) by Toth |28j is tight. If we substitute 
the current 0{k^^^n) upper bound to gtin), the time complexity is better than 
0{gk{n) logn) if fc = 0(n/ log'^ n), where c is a suitable constant. 

Another interesting question is how fast we can compute r largest maximal 
peaks for r < fc. If r = 1, Roos and Widmayer |Q gave a neat method to 
compute the maximum point in the fc-level in 0(nlogn -|- (n — fc) log^(n — fc)) 
time by using an efficient slope selection algorithm. We can compute r largest 
peaks in O(nlog'^n) -I- 0((Tn)^/^) time by combining Roos and Widmayer's 
technique and the above mentioned method for computing all the peaks. 

Finally, we investigate whether we can analogously treat some parametric 
matroids: Compute peaks in the trajectory of the largest element in the mini- 
mum weight base of a parametric matroid. In particular, the graphic matroid 
is of wide interest |l^ : Consider a weighted undirected connected graph 

G{x) with fc nodes and n edges, such that each edge has a parametric weight 
that is linear in a parameter x. Here, fc and n become the rank and size of the 
graphic matroid, respectively. Let T{x) be the minimum weight spanning tree 
of G{x) and consider the longest edge e(x) in T{x). Note that the minimum 
weight spanning tree becomes a spanning tree that minimizes the length of the 
longest edge. We call the edge e{x) the spanning bottleneck edge(SJiFi), and 
write SBE[x) and wsBsix) for e{x) and its weight, respectively. The naming 
comes from the fact that wsBsix) is the minimum value of w such that the 
subgraph of G{x) consisting of edges with weights less than or equal to w is 
connected. 
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The following problems are important in sensitivity analysis: (1). Compute 
the maximum value and the minimum value of wsBEix) for x (li I, where / is a 
given interval. (2). Compute all peaks of the trajectory y = wsbe{x). 

For example, imagine a system represented by the graph G where a link 
represented by an edge with a weight larger than a (controllable) threshold value 
becomes unreliable, and the edge weight depends on a parameter x linearly 
within an interval I. For a given subinterval J d I, we want to know the 
threshold value of the edge weights so that the graph remains connected for 
every a; G J. This can be reduced to the problem (1). Moreover, if we have 
computed all peaks in / as a preprocessing (problem (2)), we can efficiently query 
for the threshold value, provided that we have an efficient method (shown in 
Section 4) to compute wsbe{x) at endpoints of J. 

Both of problems (1) and (2) can be solved by computing the whole picture 
of the transitions of the minimum spanning trees, and the time complexity of 
the computation is 0{kn\ogn) by using an algoirthm given by Fernandez-Baca 
et al. jl^. Note that 0{hk{n)n^^^ log^^'^) time output-sensitive algorithm of 
Agarwal et al.|^ is better for some range of k, where hkin) is the number of of 
transitions of the minimum spanning tree. An 0{k^/^n) upper bound ||Tl[| and 
an n(na(fc)) lower bound are known for hk{n). 

Roos and Widmayer's method can be directly applied to the first problem. 
By using dynamic maintenance algorithms of a minimum spanning tree, 
the time complexity becomes 0{y/kn\ogn). Combined with range searching 
techniques, we improve the time complexity to 0[v?'l'^ k^l"^ + nk^^l^^. We give 
some discussion on the second problem, although theoretical improvement on 
the O(fcnlogn) time method remains open. 

2 Preliminaries 

2.1 Roos and Widmayer's algorithm 

Given a set Ti oin lines in the x-y plane, let Ck be the fc-level of the arrangement 
of 7i. Let p be a point on Ck that has the maximum y- value ymax- Without loss 
of generality, we assume that such a point is unique. For any given value a, one 
can decide whether ymax > a or not in 0{n log n) time: We sweep on the line h : 
y = a from the leftmost intersection point to the right to compute the levels of all 
intersection points on h with lines in the arrangement. If all intersection points 
are above the fc-level, a > ymax] otherwise, a < ymax- By using this decision 
method, a binary search algorithm works to compute p, and a weakly polynomial 
time algorithm with a time complexity O(nlognlogr) can be obtained, if each 
coefficient of the lines is a quotient number of integers with log F bits. Roos and 
Widmayer[^ applied an efficient slope selection method to transform the 
binary search algorithm into strongly polynomial, and gave an 0(nlog^ n) time 
algorithm. They further improved the time complexity to 0{nlogn + fclog^ k) 
for computing the minimum and 0{n logn + (n ~ k) log^{n — k)) for computing 
the maximum. 
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2.2 Range query and Matousek's point set decomposition 



We use well-known (although sophisticated) simplex range query data struc- 
tures [|| : We construct a data structure for a set S oi n points in a plane such 
that given a query halfplane H, we can compute the number of points in S 
located in H efficiently. If we spend 0{m) time for constructing the data space 
for nlogn < m < n^, the query time is 0{n/m^^^). The query can be done in 
polylogarithmic time by using 0{n/w}/'^) processors. The data structure uses 
0{m) space, although we do not discuss space complexity in this paper. More- 
over, we can query the number of points in the intersection of two (or three) 
halfplanes in the same query time if we ignore a polylogarithmic factor. We can 
also do reporting query by spending additional 0{N) time if the region contains 
points. 

Given a set H of n lines in a plane, we consider the set 'D{H) of their dual 
points: The dual point of a line y — ax — b is (a, 5). We construct a range 
searching data structure for 'D{H). Given a point p — {xo,yo), the set of dual 
points of lines below p is the set of points in 'D{Ti) located below the line 
Y = xqX ~ yo, where X and Y correspond to coordinates of the dual plane. 
Thus, we can compute the level of the point p in the arrangement of n lines by 
using half-plane range searching. Moreover, we have the highest line below p in 
the same query time. Also, we can query the number of lines which lie below 
both of a pair of query points. 

A main building block for the range query is the point set decomposition 
structure of Matousek |Q, which we also need to utilize directly (we only de- 
scribe its two-dimensional version): 

Theorem 2.1 (Matousek) Given a set S of n points in the plane, for any 
given r < n, we can subdivide S into r disjoint subsets Si {i = 1,2, ...,r) 
such that \Si\ < 2n/r satisfying the following condition: Each Si is enclosed 
in a triangle ai, and any line in the plane cuts at most cr^l"^ triangles among 
ai, ■ ■ , cTr where c is a constant independent of n and r. Such a decomposi- 
tion can be constructed in 0(n log n) time. 



3 Computing all peaks in A;-level 

We assume k < n/2 for simplicity from now on; if k > n/2, replace k hy n — k 
and exchange maximal and minimal in the statements. A key observation for 
the fc-level is that it is a subset of a union of k concave chains such that all 
concave vertices of the fc-level are vertices of these concave chains thus, a 
fc-level has at most k maximal peaks and fc — 1 minimal peaks. We want to 
compute all the local peaks in a given interval / of the x-coordinate value. 

We apply a version of parametric search paradigm |25[ |2^. However, before 
applying the parametric search, we start with a simpler "fc-branching binary 
search" method. Without loss of generality, we assume that no line in the 
arrangement is horizontal nor vertical. 
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We prepare two key-subroutines: one-shot query and peak counting: Let 
p{xo) be the point on the fc-level at the ^-coordinate value xq. Let £'f^{xo) 
(resp. i^{xo j) be the hne in the A:-level at the x-coordinate value xq — e (resp. 
Xq + e) for infinitesimally small e > 0. If xq is not an x-coordinate value of a 
vertex on the fc-level, ^^(a;o) = ^ti^o)- The above operation to compute the 
point (together with lines containing the point) on Ck at a given i-coordinate 
value is called one-shot query. One-shot query is an analogue of ray shooting Q , 
and thus the following lemma is basically well-known. The complexity q(n, m) 
given in the lemma is called one-shot query-time for the fc-level: 

Lemma 3.1 // we preprocess the lines in Ti. with 0{m) time for n\ogn < m < 
, given an x- coordinate value Xq, we can compute p{xq) , £'j^ (xq) and t^{xQ) in 
q{n,m) = 0{n/m^^^) time, and also in poly logarithmic time by using O in/ nn}^"^) 
processors. 

Proof By using the method given in the preliminary section, we can compute 
the level of any given point (zq, j/o) in polylogarithmic time by using 0(n/m^/^) 
processors. We now apply parametric searching to have the sequential time 
bound to compute the point p{xq). A parametric searching algorithm is usually 
stated as a sequential algorithm; however, it is naturally a parallel algorithm 
if we use a parallel decision algorithm and also a parallel sorting algorithm. 
We remark that we can do it easier without using the parametric searching if 
we examine the range searching method in precise; however, we omit it in this 
paper. □ 

The peak-counting is a routine to compute the number of peaks of the fc-level 
in a given interval J = [xo,a;i] of x-coordinate values efficiently. The following 
elementary lemma is essential: 

Lemma 3.2 Let /(xq) and f{xi) are numbers of positive-slope lines below or 
on the k-level at xq and xi, respectively. Then, the number of maximal peaks of 
Ck in the interval J is f{xo) — f{xi). 

Proof At-most-fc-level (the part of the arrangement below k -\- 1-level) is a 
union of k concave chains such that all concave peaks in the chains appear in 
the fc-level [^. If a concave chain among them has a peak in J, the slope of the 
chain must be changed from positive to negative. Thus, the number of maximal 
peaks within J is the difference between the numbers of positive slope lines at 
two endpoints. □ 

We remark that f{xo)~ f{xi) equals the number of positive slope lines inter- 
secting the segment between p{xf)) and p{xi) if the segment has a nonpositive 
slope. 

Lemma 3.3 For a given interval J of x- coordinate value, the number k{J) of 
peaks of Ck in J can be computed in 0{q{n,m)) time if we preprocess the lines 
with 0{m) time. Also, the number of maximal peaks can be computed in 
0{q{n, ni)) time. 



5 



Proof If we construct the dual of range search data structure for the set 
of Hues with positive slopes, the number of positive slope lines below a given 
point {xo,yo) can be computed in 0{q(n,m)) time. Hence, fixo) and f{xi) 
can be computed in 0{q{n, m)) time, and the number of maximal peaks can be 
computed by using Lemma ^ The number of minimal peaks is easily computed 
from the number of maximal peaks and slopes of the fc-level at endpoints. □ 
Now, we can apply a binary search paradigm to design a weakly-polynomial 
time algorithm. First, for the input interval /, we compute the number of 
peaks K — k{I) within the interval. The time complexity for this initialization 
is negligible, and obviously k < 2fc — 1. Next, we construct a data structure 
for the one-shot query in 0{m) time, where the choice of m will be explained 
later. Suppose that coefficients of the equations of lines are quotient numbers 
of logF bit integers. We apply K-branching binary search to find all peaks; 
At each stage of the binary search, we have at most k subintervals which has 
at least one local peak of Ck (such a subinterval is called an active interval), 
and we recursively search in active subintervals. Thus, after examining KlogF 
candidates of x-coordinate values, we can find all of the peaks. 

Proposition 3.4 All the local peaks ofXk in the interval I can be computed in 
0(^(7(71, m) log F) time. 

To make the complexity into strongly polynomial, we apply parametric 
search by using the parallel algorithm for the one-shot query given in Lemma |^ 
as its guide algorithm. We run the guide algorithm without inputting the pa- 
rameter value (in our case, an x-coordinate value), and decide the x-coordinate 
values of the peaks by using the sequential one-shot query and the counting 
algorithm of Lemma ^ as decision algorithms. The counting algorithm needs 
two values of the parameter, which corresponds to endpoints of each of intervals 
obtained by splitting active intervals in the parametric searching process. Usu- 
ally, parametric searching method solves optimization problems on monotone or 
convex functions. Here, fc-level is neither monotone nor convex, but it consists 
of K monotone fragments. Thus, while running the guide algorithm, there are 
at most K different critical parameter values to determine all the comparisons 
in the current parallel step that are necessary to proceed into the next parallel 
step. In precise, the number of different choices is the number of active intervals 
obtained breaking / by the critical parameter values found so far in the guide 
algorithm. We make a clone of the guide algorithm for each active interval. 
If the current interval is split into / active subintervals, / — 1 new clones are 
created. Naturally, we create at most k clones in our process. There is only one 
critical parameter value to determine the comparisons in a "usual" parametric 
search, and such a value can be found if we run the decision algorithm 0(log N) 
times if the guide algorithm is a parallel algorithm on TV processors. In our case, 
we run the decision algorithm O(KlogiV) times at each level. Thus, we obtain 
the following theorem: 

Theorem 3.5 All the peaks on Ck within an interval I can be computed in 
0(71 log n) + 0((Kn)^/'^) time if I has n local peaks. 
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Proof The parametric searching method gives 0{Kq{n, m)) time complexity 
apart from the 0{m) preprocessing time. We balance Kq{n,m) — 0{Kn/m^^'^) 
and m to have m — (nn)'^/^. If {kti)^^^ < n, we instead use m = nlogn. This 
gives the time complexity. □ 
Since k < 2k, we have the following: 

Corollary 3.6 All the peaks on Ck can be computed in O {n log n) + 0{{kn)'^^^) 
time. 



3.1 Computing selected peaks 

When K is large, it may be too expensive to compute all the local peaks. Suppose 
that we want to compute r largest maximal peaks in the input interval / for 
T <^ K more efficiently than computing all the peaks. This can be done by 
combining Roos and Widmayer's algorithm and the algorithm given above. 

We first run a binary search process with respect to y-coordinate values simi- 
larly to Roos and Widmayer's algorithm. At the intersection of the arrangement 
with a horizontal line y = yo, we compute intervals on the line which are below 
the fc-level (we call them semi-active intervals) in 0{n log n) time. We next com- 
pute the sum s(j/o) of numbers of local peaks in the semi-active intervals. We 
could apply our counting method of Lemma |^ for each semi-active intervals to 
compute the sum of maximal peaks in the intervals by using 0{Tq{n, m)) time. 
More simply, we can compute it in O (nlogn) time by counting the number of 
intersecting positive slope lines with the horizontal line during the sweep. In 
precise, we also need information of the arrangement at endpoints of the input 
interval / if one (or both) of them is below fc-level (we omit details). 

If s{yo) is not between r and 2t, we continue binary search on j/o- If s(yo) > 
2r, we increase yo while if s{yo) < r, we decrease yo. Thus, we can eventually 
find a value yo such that r < s{yo) < 2t. We have spent 0{n log^ n) time so far. 
Now, we search all peaks in the union of active intervals by using the method 
given in the previous section. The following theorem is easy to see: 

Theorem 3.7 We can compute r largest maximal peaks of Ck in an interval 
I in O(nlog^n) -I- 0((rn)^/'^) time. We can also compute r largest local peaks 
( including both maximal and minimal peaks ) in the same time complexity. 

Note that if we only use Roos and Widmayer's algorithm in a naive fashion 
to find r largest peaks, it would cost O(rnlog^n) time. Analogously, we can 
compute the r-smallest minimal peaks. 

Theorem 3.8 We can compute t smallest minimal peaks of Ck in an interval 
I in 0(?i log^ n) -I- ©((rn)^/"^) time. We can also compute r smallest local peaks 
( including both maximal and minimal peaks ) in the same time complexity. 

For this t smallest peak finding problem, we can modify Roos and Wid- 
mayer's method [|2^ to compute them in 0(n log^ n + kr log^ k) time. Although 
n -I- (rn)^/'^ < 2(n + kr) always holds, the time complexity is better than the 
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one in Theorem 3.1 by a polylogarithmic factor if n^/^ log"*^ n < k < n^/^ log*^ n 
for some constant c. Moreover, the algorithm does not use range search data 
structure, and hence much simpler and uses smaller data space. 

Proposition 3.9 The t smallest local peaks in an interval I can be found in 
0{n log^ n + kr log^ k) time. 

Proof First we search for a horizontal line h such that it intersects fc-level, 
and the number of peaks below it is not less than t. Such a line h can be found 
in O(nlog^n) time. Next, let v and w be the leftmost intersection and the 
rightmost intersection with the A:-level on h, and let J be the interval between 
them. If the fc-level at an endpoint of the input interval / is below h, we 
connect the point on fc-level at the a;- value of the endpoint to J with a segment 
to form a chain C with at most three segments. Let TCq be the set of lines 
in the arrangement intersecting with the chain C. The cardinality of TCq is at 
most 2fc because of Lemma |[ Finally, we find all peaks below h by using the 
r-branching binary search. Here, by using the windowing method of p^ , we 
should only take care of lines in Ho together with lines below endpoints of the 
chain C. There are at most 4fc such lines. Hence, this second step can be done 
in 0(fcr log^ fc) time. □ 



4 Bottleneck edge length in a 
parametric spanning tree 

Next, we consider the parametric spanning tree problem. Consider a connected 
graph G = {V, E) with fc nodes and n edges. Because of the connectivity, 
fc — 1 < n < fc(fc -I- l)/2. For each edge e G E, we associate a weight function 
Wei^x), which is linear on a parameter x. We assume that the arrangement 
generated by lines y — We{x) : e ^ E is simple, i.e., no three lines intersect at a 
point. We can remove this assumption by giving a symbolic perturbation. G is 
denoted by G(x) if it is considered as a weighted graph with parametric weights. 
For a given value x, we consider the minimum spanning tree T{x) of G{x). 

It is known that the number of transitions of the structure of the minimum 
spanning tree T(x)is 0{k^^^n), and all the transitions can be computed in 
0(fcn log n) time |15| . Moreover, the average edge weight in the minimum span- 
ning tree is a concave function in x, and the value of x maximizing the average 
edge weight of T{x) can be computed in O(nlogn) time |l^. 

As parametric matroid problems, the average edge weight is a counterpart of 
the average of y-values of fc lines below (or on) the fc-level. A natural counterpart 
of the fc-level itself in the minimum spanning tree is the longest (i.e. maximum 
weight) edge in the minimum spanning tree. The edge is also called the spanning 
bottleneck edge at x {SBE{x) in short), and its weight is denoted by wsbe{x). 
It is easily observed that WsBsix) is the minimum value of w such that the 
subgraph of G{x) constructed from the set of edges whose weights are less than 
or equal to w is connected. 
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It is natural and important problem in sensitivity analysis p8| to trace the 
trajectory y = wsbe{x) of the weight of SBE{x). Analogously to the fc-level, 
there are at most k maximal peaks in the trajectory y = wsbe{x). We want to 
compute peaks in the trajectory. 

4.1 One-shot query for the longest edge in MST 

We first consider efficient query for SBE{xo) at any given value of xo of the 
parameter. This query is called one-shot query for the SBE. A naive method 
is the following: First construct r(a;o) in 0{n) time, and select its longest edge. 
Instead, we use the Matousek's set partition. In the dual space, the dual points 
of n weight functions of the edges are partitioned into r subsets of size 0{n/r). 
Each subset is contained in a triangle, and 0(r^/^) triangles are cut by any 
query line. 

Accordingly, we partition the set of n edges into r subsets each has 0{n/r) 
edges. For each subset, we compute a spanning forest (irrelevant to edge 
weights) and store the connected components except singletons. Thus, each 
component has a forest with 0(min{fc, n/r}) edges. This computation can be 
done in 0{n) additional time. 

If we are given a parameter value xq, we sort 0{r) vertices of the triangles 
with respect to the inner product of them with the vector (xo, 1). We do binary 
search on this sorting list. We guess a vertex w, and consider a line t : Y = 
XqX + c which goes through v. We recognize the triangles which are below 
thus, the edges in the subsets associated with the triangles has weights less than 
c. We construct a spanning forest F of the union of forests in these subsets: 
since they have 0{rk) edges, this can be done in 0{rk) time. If the forest F is 
a spanning tree, we decide that v is too large in the sorting list, and continue 
the binary search. 

Otherwise, we consider the subsets associated with the triangles cut by (,. 
They contain 0{n/r'^/'^) points in total. We sort them with respect to the 
weights, and greedily insert them into F until we have a spanning tree. If we 
do not have a spanning tree, we decide v is too small, and continue the binary 
search. If we have a spanning tree, we decide u is a candidate, but it may be 
too large, and continue to search for the lowest vertex v satisfying the above 
condition, and return the longest edge in the tree for that v. Note that the 
spanning tree is not a minimum spanning tree in general; however, we correctly 
recognize the longest edge in a minimum spanning tree. 

This process needs 0{rk + n/r^^'^) time, and we do this process O(logr) 
times during the binary search. Thus, if we set r = (n/fc)^/^, the time com- 
plexity is 0(n^/^fc^/'^ log(n/fc)), which is slightly better than 0{n) if fc = o(n). 
By applying a hierarchical subdivision, we can further improve it: We fist start 
r — ri, and decompose the subset of size 0{rin) into r2 smaller subsets, where 
r2 = J^y^, and we further continue the refinement for = r\^J[ until r.i becomes 
below a constant. The query time becomes 

k{ri + r\^'^r2 + . . . + (rir2 . . . ri-iY^'^ri) + n/ {rir2 ■ ■ . r^)^/^. Setting n = 
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(n/fc)^/^, this enables 0{{nkY^'^) time computation for SBE{xo). Similarly 
to the case of halfspace range searching, we can combine hierarchical cutting Q 
of the arrangement to have a preprocess-query trade-off (we omit details in this 
version). Indeed, we have the following proposition: 

Proposition 4.1 // we spend 0{m) preprocessing time for n < m, we can do 
the one-shot query for SBE in 
d{n/{m/kY/'^ + k) time. 

Moreover, we will later use the following two-shot reporting query for a span- 
ning forest, which reports a spanning forest consisting of edges whose weight 
functions are below both of given two query points (a;o,yo) and {xi,yi). This 
can be done similarly to one-shot query (this is a counterpart of the simplex 
range searching if the one-shot query is a counterpart of the halfplane range 
searching). 

Proposition 4.2 // we spend 0{m) preprocessing time for n < m, we can do 
the two-shot reporting query in 
0(n/(m/fc)^/2 -I- k) time. 

4.2 Computing the maximum peak 

Let us consider the problem of computing the maximum peak in /. First, we 
straightforwardly apply Roos-Widmayer's algorithm. For a given y- value yo, we 
want to decide whether MaXxeiwsBE{x) < yo or not. We dynamically update 
the spanning forest associated with edges with weight below yo from a; = a;o to 
a; = a:i if / = [a;o,a:i]. If we find a value x G I such that the spanning forest 
becomes a spanning tree, we know MaXx^iwsBE{x) < yo- It costs 0{k^/'^) 
time to update a minimum spanning forest for insertion and deletion of edges. 
Suppose that we sweep on the line y = yo updating the minimum spanning 
forest. The line y = yo has at most n intersections with lines associated with 
weight functions, and hence the method needs 0{nk^/^) time for the decision. 
Thus, the maximum peak can be found in 0{nk^^^ logn) time. 

We try to improve the above time complexity. We subdivide the line y = 
yo into \n/s~\ intervals such that each interval contains at most s intersection 
points. For each interval li = [a:i,a:i+i], we perform the two-shot reporting 
query at (xi^yo) and {xi+i,yo)- The reported forest F is constructed from 
edges whose weight is less than yo both at a; = a;^ and x = Xi+i. If the forest 
has more than s -f 1 connected components, it is impossible that wsBE{t) < yo 
for a t e li- Otherwise, we dynamically maintain the spanning tree, where we 
contract nodes into at most s -\- 1 super nodes each of which associate with a 
connected component of the forest F. Our graph has only s edges, and hence 
the update can be done 0{y/s) time per intersection. 

Hence, total time complexity becomes 
d{ny^-\- {n/s)[n/{m/kY/'^ + k] -\- to). 

If we optimize this, we have 0{n^^'^ k^^'' -I- nfc^/'^). This is an improvement over 
0{nk'^/'^), since n<k{k + l)/2. 



10 



The minimum of wsbe{x) can be analogously computed. Hence, we have 
the following theorem: 

Theorem 4.3 For a given interval I of the parameter value x, we can compute 
both the maximum and the minimum of wsbe{x) for x G I in 0{n^^'^ k^^"^ + 
nk^/'^) time. 

We can generalize the above theorem for the truncated matroid of the graphic 
matroid to obtain the following proposition (wc omit the proof): 

Proposition 4.4 For a constant c, let wsbe-c{x) be the minimum value of 
w such that the set of edges in G{x) with weights less than or equal to w has 
at most c connected components. We can compute both the maximum and the 
minimum of wsbe-c{x) in 0{n^^'^ k^^"^ + nk^^^) time. 

4.3 Computing all local peaks for SBE 

It is desired to apply the method of computing all peaks in the fc-level to SBE 
in a parametric graph. It is known that there are at most k maximal peaks in 
y = wsbe{x). Unfortunately, we have no good method to know the number of 
peaks within an interval / = [a;o,a;i] exactly, since we do not have a property 
that is a counterpart of Lemma ^. The only known method for the authors is to 
compute the minimum spanning trees T{xq) and T(xi) explicitly, and compute 
the difference d{I) between the number of edges whose weight functions are 
positive slopes. For any disjoint set of intervals, the sum of d{I) over the intervals 
is at most k, and d{I) gives an upper bound of maximal peaks of y = wsbe{x) 
within /. However, it is often an overestimate, since d{I) gives the number of 
maximal peaks of k trajectories of weights of all edges (not only maximum one) 
in the parametric minimum spanning tree, where we only include peaks where 
an edge in the spanning tree is replaced by an edge outside the spanning tree. 

By using the above observation, we have an 0{kf{n, k)) time algorithm 
if we have an algorithm to compute T{xo) for a given xq in 0{f{n,k)) time. 
Unfortunately, we only have f{n,k) — 0{n), which leads to an 0{kn) time 
algorithm, which is inferior to a known algorithm to compute all transitions of 
the parametric minimum spanning tree. We can compute r largest peaks in 
transitions of edge weights in MST in 0{Tn) time, if we include all peaks of all 
edges in MST; however, the number of peaks appeared at the transitions of the 
longest edge among them may be much smaller than r. 

Although the above method is not attractive for the minimum spanning tree, 
the method is applicable to any parametric matroid, and hence it is useful if 
we do not have a dynamic algorithm to maintain a minimum weight base. The 
current 0{kn\ogn) time algorithm to compute the transitions of parametric 
minimum spanning tree needs 0{k^^^) time method (indeed, it can be done 
in 0{k^^^) time) to update a minimum spanning tree. Thus, the analogue of 
the algorithm needs 0{k^^^nq + kn) time to compute all the transitions of a 
parametric matroid of rank fc, where q is the time complexity to update its 
minimum weight base. 
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5 Concluding remarks 



The number of maximal peaks in a fc-level is known to be at most fcC^ (number 
of combinations choosing d elements from k elements) if we have d dimensions 
Hence, this is much smaller than the complexity of whole arrangement, 
especially if k is much smaller than n. However, to the author's knowledge, the 
problem of computing peaks in the fc-level for a higher dimensional arrangement 
is open. Although the algorithm of Roos and Widmayer can be applied to 3- 
dimensional case, it needs 0{n^) time to compute the largest peak (i.e. global 
maximum) if we naively implement it. One necessary constituent is to develop a 
counterpart of Lemma ^: Given an arrangement of n hyperplanes in the three- 
dimensional space, preprocess it, and for any given three points A, B, and C 
in the plane z = 0, decide whether the triangle ABC contains (a projection) 
of a peak in the fc-level or not efficiently. For the purpose, we probably need a 
counterpart of Lemma |^: Give a criterion of the existence of a peak from the 
information of the set of hyperplanes below fc-level at each of A, B, and C. 
In two-dimensional space, the lines are classified into positive slope lines and 
negative slope lines, while this kind of natural discrete classification of planes in 
the space does not exist. Moreover, it is difficult to find a counterpart of concave 
chain decomposition for three dimensional fc-level |^. These lacks make the 
problem difficult, although the authors think it is quite attractive. 

Another interesting problem is an an extension of parametric SBE problem 
to the c-edge-connectedness for c > 2. Here, we hope we can develop efficient 
solutions by combining geomertic methods and graph theoretical methods [p^, 
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